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Response Functions for a Granular Fluid 

James Dufty 

Department of Physics, University of Florida, Gainesville, FL 32611 USA 

The response of an isolated granular fluid to small perturbations of the hydrodynamic 
fields is considered. The corresponding linear response functions are identified in terms 
of a formal solution to the Liouville equation including the effects of the cooling reference 
state. These functions are evaluated exactly in the asymptotic long wavelength limit and 
shown to represent hydrodynamic modes. More generally, the linear granular Navier-Stokes 
equations for the response functions and related Langevin equations are obtained from an 
extension of Mori's identity. The resulting Green-Kubo expressions for transport coefficients 
are compared and contrasted with those for a molecular fluid. Next the response functions 
are described in terms of an effective dynamics in the single particle phase space. A closed 
linear kinetic equation is obtained formally in terms of a linear two particle functional. This 
closure is evaluated for two examples: a short time Markovian approximation, and a low 
density expansion on length and time scales of the mean free time and mean free path. The 
former is a generalization of the revised Enskog kinetic theory to include velocity correlations. 
The latter is an extension of the Boltzmann equation to include the effects of recollisions 
(rings) among the particles. 

§1. Introduction 

It is generally accepted that many states of driven (activated) granular matter 
admit a macroscopic description analogous to that of molecular fluids. Such states 
will be referred to as "fluidized". As with normal fluids, that continuum descrip- 
tion can be simple or complex depending on the particular fluid considered and its 
state (e.g., large or small spatial gradients; near or far from boundaries; laminar or 
turbulent; homogeneous or heterogeneous). 1 )' 2 ) In all cases, the basis for the macro- 
scopic continuum description is the exact local balance equations for the relevant 
hydrodynamic fields, and suitable approximations for the fluxes in these equations 
(constitutive equations). These approximations may be based on extrapolations 
from experimental observations, constrained by symmetries and other known theo- 
retical properties. The simplest case is that of a one component system for which the 
macroscopic fields are the density, temperature, and flow velocity. If in addition, the 
dimensionless spatial gradients and time derivatives are small then a Navier-Stokes 
hydrodynamics is expected to apply. However, the form of these equations and the 
dependence of their coefficients on grain properties and fluid state conditions are not 
known a priori. For granular gases, the hydrodynamic equations and explicit ex- 
pressions for their coefficients can be derived from the granular Boltzmann equation. 
The origins of the Boltzmann equation 3 ) and a hydrodynamic description 4 ) can be 
made quite clear in this case. 

At higher densities, still within the constraints of small gradients, the methods of 
linear response can be applied to obtain the linear Navier-Stokes equations and for- 
mally exact expressions for the transport coefficients; 5 ) the generalization to obtain 
non-linear hydrodynamics in this way has also been described. 6 ) The application of 
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linear response parallels closely its application for the same purpose to molecular flu- 
ids, 7 )' 8 ) although significant differences appear due to the inelasticity of collisions for 
granular fluids. The next two sections give a brief review of the central ideas to allow 
the reader to see these similarities and differences. For example, an new concept of 
global invariants is observed with a direct connection to long wavelength hydrody- 
namics at the microscopic level. A new formal approach to writing equations for 
the response functions is given in Section 4 from an extension of the Mori projection 
operator method. Navier-Stokes hydrodynamics and expressions for the transport 
coefficients are extracted from the long wavelength form of these equations. In the 
second half of this presentation the practical evaluation of the response functions is 
considered in the context of kinetic theory. As observed elsewhere, 9 ) a non-trivial 
generalization (e.g., retaining velocity correlations) of the phenomenological Enskog 
kinetic theory follows directly from a short-time (Markovian) approximation. Next, 
the results of a systematic low density expansion are described for leading order 
corrections to the Boltzmann limit. The relevant length scale chosen in this expan- 
sion is the mean free path, leading to the inclusion of "ring collisions" at this order. 
Some outstanding problems and opportunities for the future are mentioned in the 
summary section. 

The results described here are a personal perspective developed primarily by 
the author and his collaborators. The references reflect this bias, and apologies are 
offered to those whose excellent contributions to related topics are only recognized 
in the books and reviews cited here. 

§2. Liouville equation and hydrodynamic modes 

Attention will be restricted to the simplest realistic model of a granular fluid: 
a system of N smooth, inelastic, hard spheres of mass m, diameter a, and positive 
normal restitution coefficient a < 1. The state of the system is specified by distri- 
bution function p(r) for the N particle phase space r = {qi, ..qjv, Vi, ..vjv}- The 
evolution from a given initial state is given by the Liouville equation 

d t p(r,t) + Lp(r,t) = o. (2-i) 

The operator L generates the hard sphere Liouville dynamics for the N particles: 
linear trajectories until a given pair is in contact, followed by an inelastic change in 
that pair's velocities, and that process continued. On each collision there is a loss of 
energy so the total energy for an isolated system decreases in time. 10 ) It is usual to 
define a temperature in terms of that energy by E(t) = 3NT(t)/2. This is simply a 
definition and has no implications for an underlying thermodynamics. The decrease 
of T(t) is represented by a cooling rate ( = —dt InT. 

For an isolated system a special homogenous solution is sought in which all of 
the time dependence occurs through the energy, or temperature 

p h (r;n,T h (t),U) = (£v (t))- 3N pi > (2-2) 

where I is a characteristic length scale (e.g., mean free path) and vo(t) = y / 2T/ i (t)/m 
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is the thermal velocity. This is known as the homogeneous cooling state (HCS), and 
is parameterized by the average density n, the temperature Th(t), and the average 
system velocity U a , all of which are uniform. The r dependence of the dimensionless 
distribution p* h follows from the Liouville equation 

1 N 

Cp* h = 0, CX = LX + -( h J2Vv g -(VpX). (2-3) 

p=i 

with V Q = v„ - JJ a . In the definition of C the cooling rate Ch must be determined 
self-consistently, i.e. 



air*KCpl = => Ch = J dr*Kp* h , (2-4) 

where K is the total kinetic energy and K is its time derivative. 

For more general states a similar dimensionless distribution can be defined 

P (r,t) = (iv (t))- m P * ({t'^^t} ' n ^' s ) • (2 " 5) 

A corresponding dimensionless time s, C and Ch are given by 

r^^-c, a = -^Ch (2-6) 

« Vh{t) v h {t) 

If £ is the mean free path, then s has the interpretation of the average number of 
collisions in the interval (0,i). The Liouville equation for p* then becomes 

d s +T)p*=0. (2-7) 

The definition of C is such that its dimensionless form, C* , is independent of s. In 
this representation of the Liouville equation, it is seen that the HCS is a stationary 
solution. 

It is possible to identify a set of important eigenvalues and eigenvectors of the 
generator of dynamics C by differentiating (|2-3p with respect to its parameters 
n, T, U. Actually, it is convenient to change variables from nh and Th to rih and 
Ch (nh, Th). The result is the set of five eigenvalue equations 

T$ a = \* a $ a , (2-8) 

with the eigenvalues and eigenvectors given in terms of p* h 

K={^\ct-\ci-\ci-\ci ) j (2-9) 
^-a^'fc' ^ 2 "5a k ' '' :, - :r ' ^77 • (2 ' 10) 
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The significance of this lies in the fact that these eigenvalues are the same as those 
of the macroscopic balance equations for average number density, energy density 
(or temperature), and momentum density (or flow velocity) in the long wavelength 
limit. In this sense they represent hydro dynamic modes at the level of the Liouville 
equation. This is not so surprising since the corresponding modes in the elastic limit 
are those of the corresponding conserved quantities, all with eigenvalues 0. The 
elastic limit results are recovered since in that limit — > and pt — > equilibrium 
distribution. 

§3. Linear response 

For molecular fluids Onsager's regression hypothesis states that the dynamics 
of spontaneous fluctuations in an equilibrium fluid are the same as those for the 
approach to equilibrium for weakly perturbed non-equilibrium states. This has been 
formalized by the extensive theoretical representations of linear response, and their 
exploration by simulations and experiment. 7 ^' 8 ) An extension of Onsager's observa- 
tion to granular fluids is described here, with special attention to hydrodynamics for 
the large space and time scale non-equilibrium dynamics 

3.1. Response functions 

The objective is to describe the dynamics of the hydrodynamic fields due to 
small initial perturbations of the HCS. These fields are chosen to be the number 
density n (r, t), cooling rate £ (r, t), and flow velocity U (r, t) 

y a (r, t) O {n (r, t) , ( (r, t) , U 3 (r, t) , U 4 (r, t) , U 5 (r, t)} . (3-1) 

The temperature has been replaced by the local cooling rate, ((n,T) = c(n)T 1 / 2 , 
since the dimensionless macroscopic balance equations are diagonal in this represen- 
tation in the long wavelength limit. The deviation of these fields from their values 
in the HCS are given by 

5y a (r, t) = j oira a {T; r) (p (T, t) - p h (F, t)) , (3-2) 

where the phase functions a a (T; r) are linear combinations of those for the local 
density, energy density, and momentum density; their dimensionless forms are given 
below. The initial state is chosen to be a "local" HCS in the sense that the state is 
perturbed only through the values of these hydrodynamic fields fields (see (|2-2p ) 

«fc({v«(0)}) = 11^^'°))^ ( { qp /-- U ^°) },n(q a , )/3 N ) . (3 . 3) 

V 1 « (q a ,0) J 

This can be viewed as a state for which each small cell of the fluid is in a HCS asso- 
ciate with the local values {y a (r,t)}. It is not a solution to the Liouville equation, 
but rather a state whose average fields are specified as functions of r, t. In this sense 
it is the analogue of the local equilibrium ensemble for molecular fluids. 
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Use of this initial state in (|3-2p and expansion of pih about the uniform values 
of {y a } for the HCS to linear order gives the linear response 

5y a (r, t) = J dv'R a p (r,t; r', 0) % (r', 0) . (3-4) 

The response functions are identified as 

R aP (r, r', t) = J dra a (r, r) e - Tt ^ (r ; r') . (3-5) 

with 

(r; r) = g J(qu " r) ^(q„o) ■ (3 - 6) 

It is seen that the functions {(ftp} are just the densities associated with the hydro- 
dynamic modes of (!2-8p . J drcftp (T\ r) = <P a (r). This is analogous to the linear 
response of a molecular fluid to an initial local equilibrium state, where the pertur- 
bations are proportional to the densities associated with the conserved quantities 
(the long wavelength hydrodynamic modes). 

It is appropriate at this point to introduce the dimensionless variables of the last 
section. Also, since the reference HCS is uniform the response function depends on 
r, r' only through r — r' and a Fourier representation is convenient 

/(k*) = f drV k *- r 7(r*), /(r*) = -L £ e^*** f(k*). (3-7) 

J k* 

The dimensionless response functions are then 

R* a/3 (k*;s) = ^J dr*a* a (r*X)e- rs ^(r*;-k*) (3-8) 

with 

a* a (r*X) = £ e ik *' q ^ (v*) , a a (v*) & jl, ^ - \ + £* ■ v*,^ • v*,e 2 • v*} 

(3-9) 

It follows from the definition of i?* /3 (k*;s) at s = that {a a } and form a 

bi-orthogonal set in the sense 

i | dT*2*(r*,k*)^(r*;-k*) = <J Q/3 . (3-10) 

3.2. Hydrodynamics 

If the fields obey a closed set of hydrodynamic equations on large space and 
time scales, that dynamics must also be reflected in the dynamics of the response 
functions (Onsager's regression relationship) 

R* (k*; s) -> ^(k*)e- A * (k * )s , (3-11) 
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where A* (k*) is matrix whose eigenvalues are those for the hydrodynamic modes. 
The matrix can be given a precise representation in terms of the response functions 

A* (k*) = — lim(d s R*)R*~ 1 . (3-12) 

The limit indicated means long times and small k* . A straightforward expansion 
to order k* 2 leads to the Green-Kubo expressions for granular Navier-Stokes hy- 
drodynamics. This will not be repeated here, other than to note from (|3-6p that 

^<f>p (jT*; — 0*)| are the long wavelength hydrodynamic modes of the Liouville equa- 
tion so it follows directly that A*^ (0) = A*<5 Q/ g, 

R* af3 (0; s) = e- x >5 a p, (3-13) 

with A* given by (pU)) . 



§4. Mori identity and generalized hydrodynamic equations 

Formally exact equations for the response functions can be obtained by the 
projection operator method developed by Zwanzig and Mori. 11 ) Its generalization to 
the granular response functions here is straightforward. To simplify the notation a 
matrix representation will be used and the asterisk and dependence on k suppressed 
for the moment. Define the bi-linear functional (x, X ) by 

(x, x ) = y J drx(r) x (r), (4-i) 

where x and x are from the dual spaces of a and 4>, respectively. It is understood 
that dimensionless variables are used throughout so the asterisk is also suppressed. 
The response function is then 

R(s) = (a(s),<l>) = (a,ct>(s)), (4-2) 

with the phase space and Liouville dynamics defined by 

a (s) = e Cs a, (s) = e~ Zs (/>. (4-3) 

In this context C is the adjoint of —C. 
Define the projection operators 

Vx = (x, (/>)<!, Vi x = <f>(a, x ). (44) 

From (EHUD , (a,0) = 1 and it follows that V 2 = V, V^ 2 = T>\ and (Vx,x) = 
[x, P+x)- The response function has the representations 

R(s) = (Va(s), ( f>) = (a,vU(s)). (4-5) 



Thus the response functions depend only on the projected dynamics. Mori's identity 
is a decomposition of the time evolution of a (s) or cp (s) into contributions from this 
projected dynamics and a remainder orthogonal to it (see Appendix) 

a(s) = R(s)a+ [ drR(s - r)h (r) , h (s) = e QCQs QCa, (4-6) 
Jo 
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0(s) = 0R(s)+ f dTj(s-T)R(T), 7 (s) = -Q)e~ Q]CQU Q}C(t>, (4-7) 
Jo 

where Q = 1 — V. The equation for R(s) now follows directly from differentiating 
^2|) and use of Mori's identity 



d s R(s) + QR(s) + [ drA(s - t) R(t) = 0, (4-8) 

The matrices generating the dynamics are 

f2=(a,£(f), A(s) = -(/i(0),7(s)). (4-9) 

Finally, differentiating the identities and use of (|4-8|) leads to new equations for 
a(s) and <p(s) 

d s a(s) + f2a(s)+ I drA(s - r) a(r) = h(s) , (4-10) 
Jo 



s 



d s cf>(s)+(t>(s)f2+ dT^(r)A(s-T) =7(5). (4-11) 







The left sides of the equations represent the macroscopic or average dynamics of the 
response functions. The functions h (s) and 7 (s) on the right sides are orthogonal to 
a(s) and 4>(s), respectively, and have the interpretation of sources for "noise" . Conse- 
quently, these equations have the interpretation of generalized Langevin equations. n ) 
In this context, the second equality of (|4-9j) represents a "fluctuation-dissipation" re- 
lation. All of the results to this point are still exact. 

4.1. Hydrodynamics 

The relationship of the formally exact equation (|4-8p to hydrodynamics is most 
easily discussed in terms of the Laplace transform of the response functions 

/•oo ^ 

R(k,z) = dse~ zs R (k; s) , ReOz 0l (4-12) 
Jo 

where it is assumed that there exists some real zq such that R a a (k; s) < e~ z ° s for 
all s. Then (|4-8p gives 

R(k,z) = (z + f2(k) + A(k,z))~\ (4-13) 

where A(k, z) is the Laplace transform of A(k, s). The zeros of det (z + i?(k) + A(k, z)) 
define the spectrum of R(k, z) and the complete dynamics of R (k; s). The spectrum 
for k = is known from the exact result (|3-13p showing a set of poles at z = —X a . 
Assuming only continuity in k these poles shift at finite k to the values A a (k) de- 
termined from 

det (-Aa (k) + fl(k) + A(k, — A Q (k))) = 0. (4-14) 
If the dependence on k is analytic then the solutions to (|4-14p can be constructed as 

a power series, A Q (k) = X a +kXa + k 2 Xa +.. The eigenvalues truncated at order k 2 
define the Navier-Stokes excitations. 
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The corresponding hydro dynamic equations are those of (|3-lip 

(d s + X(k))R(k;s)=A(k), (4-15) 

with the identification of A (k) as the residues of the hydrodynamic poles and 

A(k) = ft(k) + A(k,-A(k)) (4-16) 

The explicit Navier-Stokes form is obtained by expanding the expressions for i?(k) 
and A(k, -A (k)) from ([4^]) to order k 2 

A(k) -> A + jfeflW +A; 2 (tt (2) +A (2) (z = -A)) . (4-17) 

The terms of through order k define the Euler hydrodynamics for a granular fluid, 
while the contributions of order k 2 define the Navier-Stokes order transport coeffi- 
cients. The latter are the Green-Kubo expressions. 

There are several important similarities and differences between the hydrody- 
namics for granular and molecular fluids. At the Euler level, there are the excitations 
A which lead to a long wavelength instability not present for molecular fluids. The 
hydrostatic pressure has a dependence on temperature and density determined by 
the reference HCS distribution rather than the Gibbs distribution for a molecular 
fluid. Finally, there are dissipative contributions from the dependence of the cooling 
rate on the divergence of the flow field - a new transport coefficient. 

At the Navier-Stokes level the viscous dissipation has the same form in both 
cases. However, Fourier's law is modified by an additional dependence of the heat 
flux on the gradient of the density. This can be traced to the failure of Onsager's 
reciprocal relations, or time reversal invariance, in the granular fluid. Finally, there 
are important differences in the form of the Green-Kubo expressions for transport 
coefficients. As an explicit example, the shear viscosity from (I4-17P is found to be 5 ) 

V = ~^ J dr T xy M xy + hm J ds' J drT xy e~ s ( C+ ^)r xy , (4-18) 

where the Tjj is the usual volume integrated momentum flux tensor, and the new 
flux T xy and moment Ai xy 

r xy = -$(c + ^jM xy , M *y = -\Ys(^^ y + ^,v-^)ph- (4-i9) 

There are two parts, a time integral of a flux-flux correlation function as in the 
usual expressions for molecular fluid, and a time-independent contribution. Both 
terms occur in the elastic limit, with the latter being a peculiarity of the hard 
sphere interaction. However, for the granular fluid this expression occurs even for 
continuous particle interactions and is due to the absence of time reversal invariance. 
Also in the elastic limit the flux T xy — > T xy and the correct Green-Kubo expression 
for the viscosity of a hard sphere fluid is regained. However, the granular fluid flux 
T xy ^ Tij . Instead one of the fluxes is associated with the local conserved densities 
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a a of (|3-9[) while the other flux is associated with the bi-orthogonal densities <pp of 
(|3-6p . Finally, it is noted that the flux- flux correlation function has the invariant 
subspace of the dynamics projected out, a necessary condition for the convergence 
of the time integral for large s. 



§5. Kinetic theory of response functions 



Practical evaluation of the response functions for all k, s is a formidable task for 
both molecular and granular fluids. In the former case, phenomenological approaches 
based on modeling "memory functions" in terms of a few moments has met with 
considerable success. 12 > More controlled and detailed descriptions are provided by 
kinetic theory. In this section the representation of response functions appropriate for 
such methods is given. 9 ) Subsequently, two practical approximations are described. 

The Fourier transformed response function can be written in the alternative form 

R a p(k;s) = J dra a (r,k)e- Zs (f>p(r;r = 0). (5-1) 

Then, since the local densities a a (r, k) are sums of single particle functions (see 
(!3-9j) ). integration over all degrees of freedom except those for one particle can be 
performed to get a representation in the single particle phase space 

R aP (k; s)=n j dvidqie ik ^a Q (vi) <j>$ (x i; r = 0, s) , (5-2) 

The reduced density 4$ (xi; 0, s) is a member of a set defined by 

n m ^ (si,---, x m ; 0, s) = — — [ dx m+1 ..dx N e- Zs cl> a (P; 0) . (5-3) 

(TV — my. J 

Here Xi = q^, Vj denotes a point in the six dimensional phase space for a single parti- 
cle. These functions obey a hierarchy of equations obtained by direct differentiation. 
The first of these equations describes the dynamics of the single particle function 
<f$ (zijO, s) 

f0, + ic(3 + v 1 -V vl )+vi.Vq 1 W> =n(j2 f dx 2 T(l,2)4>M{x x ,x 2 ,s), (5-4) 

where T(l, 2) is the binary collision operator describing velocity changes due to hard 
sphere interactions (its detailed form 13 )' 14 ) is not needed here). Equation (|5-4p shows 
the coupling of 4>a to the two particle function <j>a \ which in turn obeys an equation 
coupled to 4>a ; hence the structure of a hierarchy. A closed kinetic equation for (pa 

(2) 

results from this first hierarchy equation if <fi a can be expressed approximately as a 
functional of (f>$ 

<j>W(xi,x 2 ,s) ~ <P(x 1 ,x 2 ,s | </>$). (5-5) 



Use of this functional in (|5-3p gives the desired kinetic equation, and its solution 
gives the response function via (|5-2[) . 
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Generally, such a functional relationship entails limitations of space and time 
scales, density or other conditions for validity. An example is the Markovian ap- 
proximation. This represents the functional by its exact form at s = 



$(xi,X2,S = H 



(5-6) 



This is a type of mean field approximation, representing the collisional process by 
its average value in the initial state of the system. It is therefore exact at short 
times but only approximate thereafter. For the hard sphere interactions considered 
here, it describes the non-trivial binary collisions, modified by correlations with the 
environment. It neglects dynamical correlations that develop over longer times. The 
detailed form of the Markovian approximation is worked out elsewhere 9 ) and only 
the result is quoted here. Equation (|5-4p becomes the kinetic equation 



vi + vi ■ V qi + I 



0, 



(5-7) 



where the collision operator is given by 



(i) 



dx s+1 T(l,2)<P( Xl ,x 2 ,s = + 

dx 2 T(x 1 ,x 2 )g%\x 1 ,x 2 ) \fi 1 \v 1 )d>^^2,s) + fi 1) (v 2 )^\xi,s) 



I 



dq2CA(vi,qi 2 )— / dv 2 a\ (v 2 )#(i2,s) 
nh 



(5-8) 



(2) 

The function g h (x±, x 2 ) describes the pair correlations in the HCS (including veloc- 
ity correlations), and c,\(vi,qi2) is related to its functional derivative with respect 
to the local hydrodynamic fields y\{<\ 2 ) 



CA(vi,qi 2 )= / dxT^x^M/rM 



e(l)/ 



8gf^[xi,x\5y\ 

$y\(<i2) 



(5-9) 



5y=0 



The corresponding response functions are determined from the Fourier transform of 
this equation, with the result 



R a p (k; s) = nf dv iaa ( Vl ) e^"^^-^^).^ 



(i) 

lh 



\Sy=0, 



(5-10) 



where 1(h) is the Fourier representation for the collision operator /, and the initial 

condition is the functional derivative of the single particle distribution deter- 
mined from the local HCS distribution f)3-3[) . 

The Markovian approximation leads to a practical result, without a priori lim- 
itations on the density, or length scale. Thus it can be applied to both dilute and 
dense fluid conditions, from hydrodynamic length scales to those smaller than the 
grain size. In the elastic limit it reduces to the generalized Enskog approximation 
that has met with great success in that diverse context. 7 )' 12 - 1 For Granular fluids, 
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(2) 

the required input property g h (x\,x 2 ) is not known and hence most applications 
to date have used the further approximation of neglecting the velocity correlations. 
Still, the results provide an accurate and practical extension of the Boltzmann results 
to finite density gases. 

§6. Response at Low Density 

A second example of a kinetic theory description for response can be obtained 
more systematically using a dimensionless density as expansion parameter. The 
starting point is the hierarchy for the reduced distribution functions (x\ , • • • , x m , s) 
associated N particle distribution function evolving from the initial local HCS dis- 
tribution (13-31) 



n m / (m) • • ■ ,x m ;s) = (Ar _ m) , J dx m+1 ..dx N e' Cs p eh ({y a (0)})- (6-1) 

It follows that the function (/)& ^ (x\; r = 0, s) that determines the response functions 
in (|5-2[) is obtained from f^\xi,s) by functional differentiation 



( Xl ;r = 0,s) 



5f^(x 1 



(6-2) 

5y=0 



The analysis proceeds by first finding a kinetic equation for /W and then determining 
the corresponding kinetic equation for 4>a by functional differentiation. Only an 
outline of the results is provided here, with further details 15 ^ to be given elsewhere. 
The first hierarchy equation is 16 ) 

^ s + icV Vl -v 1 +v 1 -V qi ^)/«(x 1 ; S ) = J dx 2 T(l,2)fW(x 1 ,x 2 ;s). (6-3) 

A kinetic equation is obtained by finding a functional relationship f^ 2 \xi,xj;s) = 
F^(xi,Xj,s | f^). This is accomplished as an expansion in the parameter o~/£, 
where a is the grain diameter and I = 1/na 2 is the mean free path, for o~/£ << 1. 
The hierarchy equations can then be solved formally for /( m ) as a power series in 
this small parameter. The solutions can be inverted to give a representation of /( m ) 
in terms of /W and the initial conditions. The result is 



FW( Xl ,x 2 ;s I /W) = / (1) (x 1 ; S )/( 1 )(x 2 ; S ) + (~) C^(x 1 ,x 2 ;s\ /«) + .. (6-4) 
where is determined from 

(d a + £(l | +C{2 | /«)) C {2) ( Xl ,x 2 ;s | /«) = T(l,2)fW(x 1 ;s)fV(x 2 ;s). 

(6-5) 
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The single particle generator C (l | f^- 1 ') is 

c (1 1 /«) x(i) = (vi ■ v qi + icv vi • Vl ) X(l) 

-J dx 3 T(l,3) ( K f^(x 1 -s)X(x 3 ) + X(x 1 )f^(x 3 ;s)) .(6-6) 

These results determine the functional F^> (x\, x 2 ; s | •) exactly to order (a/£) 2 . 
For the special case 5y a = in (|6Tp . these equations become time independent and 
determine the one and two particle reduced distribution functions for the HCS at low 
density. Closely related equations have been studied to determine pair correlations 
in the HCS 18 ) 

To expose the content of this result, the linear equation for can be solved 
CW(xi,x 2 ;s\ /W)=W(a,0) 

+ f ds'U(s, s')T(l, 2)/W (x 1 ;s')f^(x 2 ; s% (6-7) 
Jo 

where the solution operator U(s, s') obeys the equation 

(d s +£(xx\ + c(x 2 \ U{s, s') = 0, U{s',s') = l. (6-8) 
With these results, the first hierarchy equation (|6-3|) becomes 

(^ + icVv 1 -v 1 +v 1 -V qi )/«(x 1 ; S ) = J dx 2 T(l,2)f^(x 1 ;s)f^(x 2 ;s) 

+ ( j) 2 f ds' J dx 2 T(l, 2)U(s, s')T(l, 2)/« ( Xi; S ')/M (a!2 ; a ') 

+ (l)V ^(1,2)^,0) 

To lowest order, the first term on the right side gives the granular nonlinear Boltz- 
mann kinetic equation. This approximation has been studied in some detail for 
low density granular gases. 2 - 1 ' 17 ) The second term on the right side describes "ring 
collisions" , which are dynamically correlated recollisions among pairs due to sequen- 
tial collisions with many other particles. Finally, the last term describes related 
dynamical correlations due to initial static correlations in the initial local HCS. 

The linear kinetic equation for (j>a to determine the response functions fol- 
lows from (|6-2p . functional differentiation of (|6-9p . and setting 5y a = 0. It has the 
structural form 

{d s + C h (si)) <j$ On; s) + / ds'K (s - s') 0« (xi; s') = I a On; a) , (6-10) 

Jo 

where the source I a is due to the evolution of initial pair correlations, inherited from 
the last term of (|6-2p . The operator Ch (1) is the generator for dynamics of the linear 
granular Boltzmann equation 



ft , ^2)-f^(x 1 )f'>(x 1 ) 



X2)-fH\x l )fg\x 1 ) 



(6-9) 
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C h (xi) = fvi • V qi + ^C h (3 + Vl • V Vl )) X( X1 ) 

- J dx 3 T(l,3) (f£\ Vl )X(x 3 ) + X( Xl )f£\v 3 )y (6-11) 

It can be shown that Ch (1) is the low density limit of the generator for the Markovian 
approximation (|5-7p . The last term on the left side of (|6T0p represents the correlated 
collisions generated by the linear ring operator TZ (s), and the source on the right side 
is due to the dynamical evolution of initial pair correlations. Both of these terms are 
proportional to (o~/£) 2 . Their detailed forms are lengthy and will not be given here. 
In the elastic limit, they reduce to the ring kinetic theory for response functions of 
a molecular fluid. 

The most interesting feature of the ring operator is its description of a mech- 
anism for slow dynamics not present at the level of the Boltzmann equation. It 
is due to a coupling of hydrodynamic modes that originates from the generator 
e -(C h (x 1 )+c h (x2})s j n so i u ti on to (|6-8p . It is well known that the hydrodynamic 
modes of a molecular fluid are contained in the spectrum of the linearized Boltzmann 
operator, and a similar result has been proved for the linearized granular Boltzmann 
equation. 19 ^ Thus in general the action of this operator generates the exponential 
decay of two hydrodynamic modes at all possible wavelengths, leading to slow de- 
cay (algebraic, logarithmic) of response functions for molecular fluids. Similar slow 
decay has been predicted on more phenomenological "mode coupling" methods for 
the related flux correlation functions in granular fluids. 22 ) The above kinetic equa- 
tion provides the means to investigate this mechanism in more detail and broader 
context. 

A Laplace transform of (|6-10j) and its formal solution leads to the Fourier and 
Laplace transformed response functions in the form 

R(k,js)=n y'dva(v)(z + /: / ,(k)+^(k,z)) _1 ^ 1 )(v)+7o(v,k,z)). (6-12) 

where ^^(v) is the reduced single particle function determined from the global 
invariant $ in (|2- lOj) . An important consequence of the eigenvalue equation (|2-8p for 
the Liouville equation is that these same long wavelength eigenvalues are poles of 
the resolvent operator in (|6-12p 

(£ A (0)+£(0,-A a ))#W =A«#« (643) 

Assuming analyticity in the wavevector k this assures the existence of hydrodynamic 
excitations in the spectrum of Ch (k) + TZ (k,z). Then by perturbation theory 19 ) the 
hydrodynamic modes, including effects of the ring collisions, can be explored through 
Navier-Stokes order. 

§7. Summary 

The presentation here has given an overview of linear response functions for 
excitations in the HCS or a granular fluid. The HCS is the homogeneous state of 



14 



James Dufty 



an isolated fluid, corresponding to the equilibrium state or a molecular fluid, and is 
perhaps the simplest case to consider. The response for states of more practical in- 
terest (steady shear or other driven steady states) is considerably more complex, and 
relatively little is known even in the case of molecular fluids. Here, an attempt was 
made to show the structural features that are the same for granular and molecular 
fluids, while attending to their differences in detail. 

Although the utility of these functions to describe hydrodynamic response has 
been emphasized, there application for shorter length and time scale phenomena 
should be recognized. There are many methods for practical evaluation suggested 
by experience with molecular fluids. The kinetic theory approach has been illus- 
trated here with the Markovian approximation that encompasses all length and time 
scales and a range of densities well beyond the Boltzmann limit. In addition, system- 
matic corrections to the Boltmann limit have been described for low density granular 
gases to include the novel effects of correlated many-body collisions (rings). Still to 
be explored are the consequences of this low density kinetic theory, particularly the 
consequences of hydrodynamic mode coupling in the spectrum of the ring operator. 
While this is well known for molecular fluids, there is a qualitative difference in 
the hydrodynamic modes for granular gases. Some of the modes are unstable (in 
dimensionless form) and there is the potential for a qualitative difference in the hy- 
drodynamic dispersion relations determined from R(k, z). This occurs for molecular 
fluids in two dimensions, due to mode coupling effects, and may be more signifi- 
cant for granular fluids in all dimensions due to the unstable modes. The detailed 
evaluation of (|6T2p will clarify this possibility. 
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Appendix A 

Mori 's identity 

The projection operators of (!4-4p are Vx = (x,i/i)a, V^x = 4>(a,x)- Let 
Q = 1 — V and = 1 — denote the corresponding orthogonal projections. Then 
a (s) can be decomposed as 

a (s) = (P+Q)a (s) = R(s)a + a ± (s), (A-l) 

where a±(s) = Qa(s)Q. Next operator with Q on the equation of motion (d s — C) a (s) = 
0, to get 

(d a - QCQ) a±(t) = QCVa(t) = R(s)QCa. (A-2) 
The solution is s 

a±(t)= I dTR(s-T)e Cr QCa, (A-3) 
J o 
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and its substitution into (|A-lj) gives Mori's identity (I4-6|) . A similar analysis for tp (s) 
gives the identity (|4-7p . 

Finally, differentiating the response function and using Mori's identity gives 

d s R(s) = R{s) (Ca, (f>) + / dri?(s - r) (Ch (r) , 0) 

Jo 

= - (a,~C(j)) R(s) - [ dr (a,Cj(s - t)) R(t) (A-4) 
J o 

The first equality comes from using the form (a (s) , <^>) while the second results from 
the form (a, 4>(s)). 
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